Evanescent Gain in "Trapped Rainbow" Negative Refractive Index Heterostructures 
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We theoretically and numerically analyze a five-layer "trapped rainbow" waveguide made of a 
passive negative refractive index (NRI) core layer and gain strips in the cladding. Analytic transfer- 
matrix calculations and full-wave time-domain simulations are deployed to calculate, both in the 
frequency- and in the time-domain, the losses or gain experienced by complex-wavevector and 
complex- frequency modes. We find an excellent agreement between five distinct sets of results, all 
showing that the use of evanescent pumping (gain) can compensate the losses in the NRI slow-light 
regime. 



The ability to stop and store optical pulses could usher 
in a range of fundamentally new and revolutionary ap- 
plications [ , ], but the challenges encountered in our 
efforts to stop light are formidable. This highly unusual 
state for an electromagnetic wave refers to the situa- 
tion where the lightwave completely stops despite the 
absence of any dielectric or other barriers in the direc- 
tions where it could propagate. Unlike the confinement 
or trapping of light inside a dielectric cavity, where one 
cannot spatially separate optical bits of information, a 
stopped-light structure would allow for sequentially stop- 
ping and storing spatially separated optical bits, thereby 
potentially leading to all-optical memories [1, ]. Unfor- 
tunately, stopping of light cannot be obtained with (non- 
switchable) periodic structures, such as photonic crystals 
or coupled-resonator optical waveguides, owing to their 
extreme sensitivity to disorder, which invariably destroys 
the zero-group- velocity (ZGV) point (s) [2]. With atomic 
electromagnetically induced transparency one coherently 
imprints the shape of an optical pulse into an electronic 
spin excitation, i.e. light is "stored" but not stopped be- 
cause at the ZGV point all photons are converted into 
atomic spins and light is completely extinguished [ ]. 

A method that could allow for true stopping of light 
in solid-state structures and at ambient conditions was 
first suggested in 2007 [4]. Instead of relying on peri- 
odic back-reflections or on resonances, the deceleration 
of light in this method is based on the use of nega- 
tive effective electromagnetic parameters in metamaterial 
waveguides, which cause negative Goos-Hanchen phase- 
shift steps (i.e., deceleration) in the propagation of a light 
ray. The negativity in the real parts of the bulk constitu- 
tive parameters in metamaterials originates from the re- 
sponse of deep-subwavelength elements or layers [ ], and 
can therefore be insensitive to disorder [6]. Hence, it is 
anticipated that in metamaterial waveguides light can be 
coherently decelerated and stopped even in the presence 
of disorder and surface roughness. Indeed, time-domain 
simulations [7] reveal that stopping of light inside such, so 
called, "trapped rainbow" structures, is achievable despite 
deviations from the perfect geometry (numerical "rough- 
ness" at the interfaces) [ ]. Moreover, a recent experi- 




Figure 1. Illustration of the NRI slow-light heterostructure 
considered in the analyses. Also shown is a characteristic 
snapshot (from the FDTD simulations) of the propagation of 
the excited (dominant) slow- light mode. See also [12]. 



mental work has demonstrated trapped rainbow stopping 
of light in a tapered plasmonic system in the quasi-static 
regime [9]. 

A central task in this method of stopping and stor- 
ing light is to study whether the losses associated with 
the use of negative refractive index (NRI) metamaterials 
can be overcome. Although a series of recent works [10] 
have shown that losses in active "fishnet" metamaterials 
can be compensated, the considered structures were very 
thin in the longitudinal direction, essentially being two- 
dimensional and spatially dispersive. Furthermore, the 
dispersion relations and restrictions obeyed by light in 
fishnet metamaterials are completely different from those 
in trapped rainbow NRI waveguides. As a result, it is 
not clear until now whether in such waveguides losses 
can, even in principle, be overcome in a slow-light regime 
where the effective refractive index experienced by light 
is negative - even when use is made of gain media [11]. 

In this Letter, on the basis of analytical calculations 
and rigorous numerical simulations, we show how the 
incorporation of thin layers of a gain medium in a pas- 
sive trapped rainbow heterostructure can compensate the 
losses while simultaneously preserving the negative effec- 
tive refractive index of the guided slow light. The specific 
structure considered is illustrated in Fig. 1. It consists 
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of a NRI core layer bounded symmetrically by two thin 
gain layers that evanescently feed the supported guided 
modes [ ]. The whole structure is embedded in air. 

When gain or losses are present in one or more of 
the layers of a heterostructure the supported guided 
modes become complex. Apart from the usual complex- 
wavevector (real-frequency) solutions to the character- 
istic equation, one may then also retrieve complex- 
frequency (real-wavevector) solutions [14]. For instance, 
it is well-known from prism-coupling to uniform metallic 
films that fixing the incident light frequency and sweep- 
ing the angle of incidence results in the excitation of 
complex-A; surface plasmon polaritons (SPPs) exhibiting 
back-bending in the co-k dispersion diagram. By contrast, 
keeping the incident angle constant and varying the fre- 
quency of the incident light gives rise to the excitation 
of complex-co SPPs with distinct reflectivity dips and no 
back-bending [ ]. The imaginary part of a complex- co 
solution relates to the temporal losses experienced by a 
light pulse [14]. 

A numerical framework that is well-suited for the 
study of these modes is the finite- difference time-domain 
(FDTD) method [ ]. With this method one can accu- 
rately launch the desired negative phase velocity (back- 
ward) mode into the waveguide and directly investigate 
its temporal losses [ ]. In our simulations we deploy a 
modified total-field/scattered- field formulation [15], with 
the excitation plane oriented perpendicularly to the cen- 
tral axis of the heterostructure of Fig. 1, and the ampli- 
tudes of the H z - and E y - field components along the plane 
being set to match the transverse profile of the backward 
TM2 mode. The central frequency of the injected pulse 
is fixed to 400 THz (Ao = 750 nm), the side of the square 
FDTD cell has a length of Ax = A /200 = 3.75 nm and 
the Courant value is set to 0.7. The width of the core 
layer is w c = 0.35Ao = 262.5 nm, while the width of the 
gain layers in the cladding is w g = 0.25Ao = 187.5 nm 
whenever they are incorporated in the heterostructure. 
We model the passive NRI of the core layer using a 
broadband Drude response [ , Chap. 1.3]: ud(oj) = 
1 - uj 2 p /(uj 2 + iuT D ), with uj p = 2tt x 893.8 x 10 12 rad/s 
and Td = 0.27 x 10 12 rad/s. The frequency response 
of the permittivity of the gain layer obeys a Lorentzian 
dispersion: (uj) = e^-^ Aeuo\ / (uj\ — ilV^u — uj 2 ) , with 
= 1.001, Ae = - 0.0053, uj l = 2tt x 370 x 10 12 rad/s, 
and = 10 14 rad/s, resulting in a line-shape that is 
similar to that produced by, e.g., an electronic transition 
in a quantum dot [18]. 

When the losses and gain are relatively small, as those 
used in the considered NRI heterostructure, the imagi- 
nary part of the complex- co solution is proportional to 
the imaginary part of the complex-A: solution by a factor 
that is close-to-equal to the group velocity [19]. Hence, 
in this case, one has an additional opportunity (see Fig. 4 
later on) to check the accuracy of the obtained numeri- 
cal results by examing whether such a proportionality is 
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Figure 2. Snapshots of slow-light pulse propagation along 
the central axis of the considered waveguides for: (a) case I 
(neither loss nor gain); (b) case II (loss but no gain); (c) case 
III (both, loss and gain); and (d) case IV (gain but no loss). 
The lighter the color, the later in time the snapshot is taken. 



fulfilled. Following the standard theory of active optical 
waveguides [ ], we assume that the saturation inten- 
sity for the gain medium is sufficiently large, leading to 
a correspondingly large value of the critical gain-length 
product beyond which gain depletion owing to amplified 
spontaneous emission (ASE) can become significant [21]. 
Operation sufficiently below this limit implies that we 
are in the linear regime where no gain depletion occurs, 
and that the effect of ASE on the signal gain may be dis- 
regarded, in accordance with analogous studies of active 
optical structures [20, 22, 23]. 

To analyically confirm the accuracy of the numerical 
results we have developed a frequency-domain transfer- 
matrix method (TMM) capable of identifying both the 
complex- co and complex-A: modes of the multilayer NRI 
heterostructure. The method derives the dispersion 
equation and uses the argument principle method to lo- 
cate and isolate its zeros either on the complex- or the 
complex-A: plane [ ]. Suitable conformal mappings are 
deployed in both cases to unfold the four-sheeted Rie- 
mann surfaces associated with the characteristic equa- 
tions. Upon isolation of a zero, the Newton-Raphson 
method is used to pinpoint and track its location on the 
complex plane. The technical details of the overall pro- 
cedure will be presented elsewhere. 

In our analyses we consider the following four cases: (I) 
neither loss in the NRI core layer {Yd = 0) nor gain in 
the cladding layers (the cladding is only air), (II) loss in 
the NRI core layer but no gain in the cladding, (III) both 
loss in the NRI region and gain in the cladding strips, and 
(IV) the NRI core layer is modeled as being lossless and 
gain is used in the two cladding layers. 

In order to validate the causal dynamics of the injected 
pulse with the FDTD method we created an animation 
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Figure 3. Comparison between FDTD (symbols) and TMM 
(lines) calculations of the absorption coefficient a (spatial 
losses) versus frequency for the TM2 mode in case I (red), 
case II (blue), case III (green) and case IV (orange). The 
inset depicts the frequency dispersion of Re{n e ff} in all four 



of its propagation along the considered heterostructure 
[12]. From there one may directly see that first, a sin- 
gle slow-light guided mode is excited and second, the 
mode experiences an effective refractive index n e ff with 
a negative real part, having antiparallel phase and group 
velocities. 

Sucessive snapshots of the Gaussian pulse propagating 
down the NRI waveguide for cases I-IV are depicted in 
Fig. 2. We see that for case I (neither loss nor gain) 
the amplitude of the guided pulse decays with distance 
[Fig. 2(a)]. This reduction is not a result of the pulse los- 
ing energy but arises entirely owing to group- velocity dis- 
persion, which causes the guided slow pulse to broaden, 
thereby leading to a gradual decay in amplitude. Fig- 
ure 3 (red solid line and symbols) confirms the fact 
that energy is conserved, as the absorption coefficient 
a = 2ulm{n e ff}/c (spatial losses) is zero throughout the 
frequency spectrum of the pulse. The complex n e ff shown 
in Fig. 3 for all considered cases is extracted by recording 
the amplitude of the pulse at two fixed points along the 
central axis of the waveguide over time, and then divid- 
ing the Fourier transforms of the two time series [15, 16]. 
When dissipative loss is introduced into the NRI core 
(case II) the amplitude of the pulse decreases even fur- 
ther compared to case I [Fig. 2(b)]. Figure 3 (blue solid 
line and symbols) shows that now Im{n e ff} > for all 
frequencies, as expected. 

By introducing gain in the cladding layers (case III) 
we evanescent ly pump the pulse [13] and allow for the 
compensation of its propagation losses [compare Fig. 2(c) 
with Fig. 2(a)]. Indeed, Fig. 3 (green solid line and 
squares) shows that at approximately 400 THz (central 
frequency of the pulse) the imaginary part of the ef- 
fective refractive index becomes zero, while for smaller 



frequencies Im{n eJ j} assumes negative values (amplifica- 
tion). Thus, in case III there is a continuous range of fre- 
quencies (/ < 400 THz) where we simultaneously have 
Ke{n e jj} < (inset in Fig. 3) and lm{n e ff} < (green 
line for / < 400 THz in Fig. 3). We note that the opto- 
geometric parameters of the heterostructure have been 
chosen such that a light-pulse experiences almost the 
same frequency dispersion for all cases, with differences 
between the various Re{n ej ^} (cases I-IV) being indis- 
cernible at the linear scale of the inset in Fig. 3. For 
all cases presented we find that the parameters retrieved 
from the FDTD simulations (symbols) are in excellent 
agreement with those calculated using the TMM (lines). 

To further confirm that light amplification is in prin- 
ciple possible in the negative index slow-light regime we 
"switch off" the losses, while maintaining the gain in the 
cladding strips. Figure 2(d) shows that in this case (IV) 
the negative-phase- velocity slow pulse is amplified while 
propagating along the waveguide. In particular, it is seen 
that the pulse amplitude at around x = 40 ^m exceeds its 
initial amplitude despite the fact that the pulse has been 
broadened due to group-velocity dispersion. This con- 
clusion is further confirmed by finding that Im{n eJ ^} < 
throughout the spectrum of the Gaussian pulse, as shown 
in Fig. 3 (orange solid line and symbols). 

Next, we examine how the spatial and temporal losses 
(or gain) experienced by both the central frequency of 
the pulse and the pulse as a whole vary with core thick- 
ness (Fig. 4). The complex- co solutions can be calculated 
with the FDTD method by recording the spatial varia- 
tion of the field amplitude along the central axis of the 
heterostructure at two different time points, and then 
dividing the spatial Fourier transforms of the two longi- 
tudinal spatial profiles. The rate of energy change for the 
whole wavepacket (total loss or gain) is calculated using 
the discrete Poynting's theorem integrated over a spatial 
region sufficiently wide to contain the pulse. 

Figure 4 shows that for core thicknesses above 262 nm 
the central frequency of the pulse experiences loss. For 
smaller thicknesses, for which the amplitude of the field 
increases inside the gain region, we find that the gain 
supplied by the cladding strips overcompensates the loss 
induced by the core layer. At a core thickness of 262 nm 
the central frequency experiences neither gain nor loss, 
while the wavepacket as a whole experiences gain (inset 
in Fig. 4). In all cases we have verified that Re{n ej ^} < 
(data not shown here). 

Overall, we find excellent agreement and consis- 
tency between five distinct sets of results: the spatial 
losses/gain (multiplied by the group velocity [19]) for 
the central frequency as calculated by the FDTD (green 
dots) and the TMM (green dashed line), the temporal 
losses/gain for the central frequency as calculated by the 
FDTD (blue squares) and TMM (blue dashed line), and 
the temporal losses of the whole wavepacket as calcu- 
lated by the FDTD method (purple symbols in the in- 
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Figure 4. (color online) Comparison between FDTD 
(symbols) and TMM (lines) calculations of the temporal 
losses/gain and group velocity of the complex- co and complex- 
k solutions with varying core thickness (case III). Shown are 
the group velocity (y g ) of the complex- co solutions (black), the 
group velocity of the complex- A; solutions (red dashed line), 
the imaginary part of the complex- co solutions (blue) and the 
imaginary part of the complex- A; solutions multiplied by v g 
(green). The inset shows the rate of energy loss (or gain) 
for the whole wavepacket (purple symbols) with varying core 
thickness as calculated by the discrete Poynting's theorem 
within the FDTD method. 

set of Fig. 4). This fact provides further evidence that 
loss compensation is in principle possible in the slow-light 
NRI regime. 

Finally, we note that for core thicknesses smaller than 
around 140 nm the group velocity of the complex- A; 
mode characteristically differs from that of the complex - 
co mode (red dashed and black solid lines in Fig. 4). As 
with the case of SPPs in plasmonic films [14], the group 
velocity of the complex- A: solutions exhibits a "back- 
bending", never becoming zero, while that associated 
with the complex-co solutions may reduce to zero even 
in the presence of excessive gain (or losses). 

In conclusion, by studying in the time- and frequency- 
domain the complex- co and complex- A: modes in NRI 
slow-light waveguides with gain in the cladding region, 
we have shown that it is possible to compensate the dis- 
sipative optical losses. This geometry allows for lossless 
or amplified slow-light propagation in the regime where 
the real part of the effective refractive index experienced 
by the guided modes remains negative. We believe that 
this work could aid the realization of lossless metamate- 
rial waveguides to be used in a wealth of photonic and 
quantum optics applications. 
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